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By reversing the usual direction of application, a common procedure for solving integral equations 
numerically is used to obtain the asymptotic P-condition numbers of two well-known test matrices. 
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Todd fl] 1 has recently suggested the matrix An 
= (\i~ J\)-> 1 ^ hj^ n , as an addition to the well-known 
collection of test matrices [2]. The P-condition num- 
bers of the matrices are used as a measure of their 
difficulty for numerical purposes. Where the condi- 
tion numbers are not explicitly known, the asymptotic 
behavior in n is given. Lehmer's matrix A 7 



aij 



-ilj 
-jli 



]<i 



is exceptional in that the correct order in n is not 
known. A simple idea will allow us to obtain the 
asymptotic condition number of An and the correct 
order for A 7 . 

Hilbert's first method for integral equations [3] ap- 
proximates the eigenvalues of the kernel K(x, y), 
^ x, y ^ 1, by those of the matrix 



I 



i J 



K 

n \ \n n 



If K is bounded and Riemann integrable, then the eigen- 
values of the matrices tend to those of the integral 
equation as n tends to °c. We reverse this procedure. 
We wish to estimate the behavior of the eigenvalues 
of a set of matrices as n tends to °°. If we can regard 
them as arising from the application of Hilbert's first 
method to a fixed kernel, then we may hope for an 
asymptotic result. 

To estimate the largest eigenvalue of A 17 , let us form 



ii(\i-j\)=l 
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We regard them as arising from the approximation 
of the kernel 

K(x, y)=\x — y\. 

A simple computation gives the largest eigenvalue of 
Kas 

x = h z- 2 

1 o 

where z {) is the unique real root of coth z = z. This 
equation has been studied [4] and z () is approximately 
z {) = 1.9967864. We conclude that the largest eigen- 
value of A l7 is asymptotically 

k 



ho 2 n 2 . 



It is easy to estimate the accuracy of the approxima- 
tion using the bounds of [5] but we shall not go into 
this. 

We need the reciprocal of the eigenvalue of A 17 
smallest in modulus. It is obtained as the largest 
eigenvalue of the inverse matrix. 

Let us introduce the matrices 



M 



2A~ 



N = 




2,-1 



1, 1 



and Q where f q n = q ln = q nl = 



q >»=^- 



<7 y = otherwise. 
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It is readily verified that 

AJ=h(Q-N) 
or 

M+Q = N. 

Q is obviously a positive semidefinite matrix. By a 
well-known inequality [6, p. 73], 

\{N) = KW + Q)^\{M). 

Rutherford [7] has shown that the eigenvalues of N are 
irr 



4 cos 2 



2n 



i = l, 2, . . ., n. 



and the eigenvector associated with \ X {N) is 



%.— sin sin 

3 n n 



7=1, . . ., n. 



The Rayleigh quotient furnishes a lower bound to 
Xi(M). A convenient approximate eigenvector is that 
associated with \\{N). For this vector 



(x, x) (n— 1)0*;, x) 

A little algebraic manipulation gives 

77 



sin' 



\ R (M) = UN)- 



Thus 



and we see 



From this 



and finally 



n(n-l)\ x {N) 

\ R (M) ^ \i(M) ^ \j(A0 
Xi(Af) = Ai(AO + 0(rc- 4 ). 



: 2 cos 2 ^--h0(n- 4 ) 



M^i 



PWit)- 



'r cos" 
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The same methods applied to the kernel 

K(x, y) — x\y x ^ y 
y/# y < ^ 

show that for ^7, \i = 0(l). The inverse is explicitly 
4* 3 



4i 2 — 1 



2*i-l 



i(i+l) 
2 1 -h 1 



i =7", i < 7i 



i=j=n 



-y|>i- 



Gerschgorin's theorem immediately gives 



= 0(n). 



Todd and Newman comment that P(A 7 ) ^ ft. From 
this we conclude that P(A 7 ) ~ en where c is a constant. 
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